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We study the cooling performance of optical-feedback controllers for open optical and mechanical 
resonators in the Linear Quadratic Gaussian setting of stochastic control theory. We utilize analysis 
and numerical optimization of closed-loop models based on quantum stochastic differential equations 
to show that coherent control schemes, where we embed the resonator in an interferometer to achieve 
all-optical feedback, can outperform optimal measurement-based feedback control schemes in the 
quantum regime of low steady-state excitation number. These performance gains are attributed to 
the coherent controller's ability to simultaneously process both quadratures of an optical probe field 
without measurement or loss of fidelity, and may guide the design of coherent feedback schemes for 
more general problems of robust nonlinear and robust control. 



As present-day engineering relies broadly and implic- 
itly on real-time feedback control methodology p], it is 
difficult to imagine our nascent explorations of quantum 
engineering advancing to technological relevance with- 
out rigorous extensions of core control theory to incorpo- 
rate novel features of quantum dynamics, stochastics and 
measurement. While significant progress has been made 
recently in terms of analyzing quantum feedback sys- 
tems [2HZ] and in experimental demonstrations of quan- 
tum feedback control [5HT7]. we still have a relatively 
limited understanding of systematic approaches to quan- 
tum control design and of the qualitative role of quantum 
coherence and entanglement between the plant and con- 
troller in a feedback loop. 

Within the elementary context of linear open quantum 
systems, James, Nurdin and Petersen [TH1 US] have uti- 
lized interconnection models based quantum stochastic 
differential equations (QSDEs) [20 23 to develop quan- 
tum generalizations of the traditional paradigms of rl°° 
and Linear Quadratic Gaussian (LQG) optimal control. 
While some of the most exciting potential applications 
of quantum feedback control involve nonlinear dynam- 
ics and/or non-Gaussian noises [2~ilJ2 7j . the linear set- 
ting is an essential starting point for rigorous study and 
presents crucial advantages in terms of analytic and com- 
putational tractability. 

Here we focus on a theoretical investigation of steady- 
state cooling of open quantum oscillators such as opti- 
cal and optomechanical resonators subject to stationary 
heating, damping, and optical probing and feedback. We 
work within an LQG framework as in the recent paper 
of Nurdin, James and Petersen [19] and utilize numeri- 
cal optimization together with fundamental analytic re- 
sults PQ bounding the best possible LQG performance of 
measurement-based feedback control schemes to estab- 
lish and to interpret quantitative advantages of coherent 
feedback for cooling-type performance metrics in certain 
parameter regimes. 
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Following recent convention, as in [TH [18l [19] , we will 
here refer to measurement-based controllers as "classical" 
controllers and to coherent feedback controllers as "quan- 
tum" controllers. This terminology reflects the general 
distinction that the signal processing required to deter- 
mine LQG-optimal control actions from a real-time mea- 
surement signal can be implemented by a classical electric 
circuit, while all of the hardware in a coherent feedback 
loop must be physically describable using quantum me- 
chanics (typically with weak damping). 



I. LINEAR SYSTEMS 

Quantum harmonic oscillators can be modeled as cas- 
cadable open quantum systems using the SLH framework 
[2"81 |2"9"] and the associated QSDEs. In the SLH frame- 
work, any open quantum system may be described as a 
triple: 

G=(S,L,H) (1) 

where S is a scattering matrix, L is a coupling vector and 
H is the Hamiltonian operator for the system's internal 
degrees of freedom. For a linear system with an internal 
state x, Sij is independent of the internal state, Li — 
AiX + Xi is at most linear, and H = ^x T Rx + r T x is at 
most quadratic. 

Armed with an SLH representation the most efficient 
way to simulate a linear quantum system is to solve the 
QSDEs, which represent coupled Heisenberg equations of 
motion for system operators and input-output quantum 
stochastic processes. Following the work of James, Nur- 
din and Petersen [18] we write the QSDEs for a linear 
system in the state-space form, 

dx(t) = [Ax(t) + a]dt + B da{t) 

da(t) = [Cx(t) +c]dt + D da(t) (2) 

Here x(t) gives the plant's internal variables; this is a 
Hermitian, operator-valued vector. A, B, C and D are 
real matrices; a and c are real vectors. The processes 
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da(t) and da(t) are quantum stochastic processes for the 
inputs and outputs, respectively. For convenience, we 
make them Hermitian as well; for a given port, one has 
da t = {dA l + dA\,(dAi - dA\)/i), where dA(t) is the 
quantum Wiener process [3J2 [3Tj following the ltd rule 
dAi dA\ = Sijdt. 

Defining = [xi,Xj]/2i as the commutator matrix, 
the ABCD parameters of ^ can be related to the SLH 
parameters as follows: 
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(Here 5, A, and A are real matrices which can be eas- 
ily constructed from S 1 , A and A, which are in general 
complex. J is a canonical antisymmetric matrix of the 
appropriate size. See Appendix |A}) 

To measure the performance of a given controller we 
need to define a cost function. For example, to minimize 
the plant's response to a noisy input one could minimize 
the steady-state expectation value of the excitation num- 
ber uva\. With (classical) state feedback and in the ab- 
sence of exogenous noise such a quadratic cost function 
would result in a Linear Quadratic Regulator (LQR) op- 
timal control problem [I], but in our optical feedback 
scenario with Gaussian input fields (vacuum or thermal 
noise) this becomes a quantum LQG problem [Tl [TO]. 

It is straightforward to concatenate and cascade linear 
systems once we have the ABCD models. We have writ- 
ten software in Mathematica to compute the ABCD ma- 
trices for an arbitrary linear quantum system |46j . This 
borrows many elements from the Modelica quantum cir- 
cuit toolkit of Sarma et al. 32 , and is similar to the 
QHDL framework of Tezak et al. [33]. The code com- 
putes the LQR cost function as a function of the plant 
and controller properties, and it would not be difficult 
to extend it to more general cost functions. Thanks to 
the linearity of our system, simulation is very fast: the 
complexity is polynomial in the size of the circuit, not 
exponential as is usually the case for quantum simula- 
tions, and for a simple system, it computes the LQR in 
well under 50 microseconds. 

Given a particular plant, the code is fast enough to 
perform a multivariate Newton-Raphson optimization 
scheme to find the (locally) optimal controller parame- 
ters. This is possible regardless of whether the controller 
has any particular structure - if the controller's struc- 
ture is left arbitrary, the code can simply optimize with 
respect to the controller's ABCD matrices, subject to the 
physical realizability conditions 



AQ + OA + BJB = 
QC T + BJD T = 
DJD T = J 



Control Objective: 
J = x 2 + p 2 !< = ° 01 

= 4a*a + 2 




(Noise, k„) k 3 = 0.01 . . < 0.01 






..T 




■ 








Homodyne Heterodyne 



OPO OPO, Squeezed I/O 



(4) 



FIG. 1: Optical cavity plant system with five possible classical 
and coherent feedback controllers 



that arise from the fact that time evolution should pre- 
serve the commutation relations between system and in- 
put/output fields [HI [19]. Optimizing with respect to 
an "arbitrary" controller takes longer because there are 
more free parameters, but the code is fast enough for 
each Newton step to take no more than 1.5 milliseconds 
on a standard laptop. 

We note that the classical steady-state LQG problem is 
a convex problem, and the optimal steady-state controller 
parameters can be derived via solution of algebraic Ric- 
cati equations pQ. In the quantum case, no such closed- 
from solutions are known and the realizability constraints 
Q make the landscape for numerical optimization non- 
convex |19) . Hence while we can be sure about the clas- 
sical optimality of measurement-based controllers for the 
oscillator cooling scenarios we consider, the coherent con- 
trollers we find via numerical optimization are merely lo- 
cal minima and can only be considered as candidates for 
quantum optimality. 



II. CONTROL OF AN OPTICAL CAVITY 

As a simple example of a quantum "plant" system, 
consider an optical cavity with a noisy input, Fig. [T] In 
the controller's absence, the cavity is driven by two vac- 
uum inputs (mirrors k\ and A;2, & n( i onc thermal input 
(mirror k%). Any noise process that is much broader spec- 
trally than the cavity linewidth can be approximated as 
a "white noise" thermal input. Without such noise, the 
cavity's internal mode decays quickly to the ground state. 
The objective in this control problem is to minimize the 
effect of the noise on the cavity's internal state - in other 
words, to minimize the photon number (cfia) of the cav- 
ity. We accomplish this by sending output 1 through a 
control circuit and feeding the result back into input 2. 
This is an LQG feedback control problem. 

Five possible controllers are shown in Figure [T] The 
classical controllers work by measuring a quadrature from 
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FIG. 2: Experimental realization of an 0P0 with a cavity 
and a nonlinear crystal. 



the cavity's output (or in the heterodyne case, splitting 
the beam and measuring two different quadratures), and 
applying a feedback signal based on this measurement 
and the controller's internal state. The "trivial con- 
troller" works by feeding the output directly back into 
mirror 2 of the plant, perhaps with a phase shift. If the 
light reflecting off of mirror 2 is in phase with the light 
leaking out of the mirror, the light lost through both mir- 
rors interferes constructively, reducing the control objec- 
tive (at a) (see also [2"rJ]1. 

The remaining two controllers shown in the figure are 
coherent controllers with memory. Unlike the trivial con- 
troller, the control signal is a function not only of the in- 
put field, but also the input's history. But unlike the clas- 
sical controllers, the input field is not measured; instead, 
it is coherently processed and the result is fed back into 
the plant cavity. These designs use an optical parametric 
oscillator (OPO, as in Fig. [2]) to squeeze the optical field. 

The OPO will have the following SLH model: 
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(5) 



Here, k,\ is related to the input/output mirror reflectance, 
K2 to other losses; A is the cavity detuning, and e is a 
complex number, whose amplitude encodes the strength 
of the pump field and the nonlinear medium's x^ 2 \ and 
whose phase encodes the pump field's phase [54] . 

The plant system, an optical cavity with a noisy input, 
can be modeled as an open quantum system with three 
couplings, one for each mirror. The SLH model for this 
system is: 



5 = 1 
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(6) 

We also need to find the covariance matrix Fij for the 
noisy inputs ddi, defined by | (daidaj + dajdai) = F^dt. 
Recall that, for vacuum inputs, the fields dA and dA' 
satisfy the Ito relations dAdA = dA^dA^ = dA^dA = 0, 
dAdA^ = dt [22 E], leading to the Ito tables: 



dX/dY 


dA dA* 
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For a non-vacuum, thermal input, the field dA has ad- 
ditional (unsqueezed) noise, so dA^dA = k n dt for some 



Linear 


SLH 


Model . 


n sys 2 , rip or t 


— 3 



















1 
















o 


/ ki 
V i 







1 

























1 















S = 












L = 

















1 


























4^ 
















1 






















1 

























, A 
2 




Inputs : dai 
Outputs: dbi 



da 2 
db 2 



da3 
db 3 



FIG. 3: Output from our Mathematica package, describing 
the plant system. 



noise strength k n > 0, and the rest of the relations are 
adjusted accordingly, leading to the following Ito tables: 



dX/dY 


dA dA^ 


dx/dy 


da x da p 


dA 


(l + k n )dt 


da x 


(l + 2k n )dt idt 


dA* 


k n dt 


da p 


-idt (l + 2k n )dt 



In the present system, inputs 1 and 2 are vacuum, and 
3 is thermal noise. This gives the following covariance 
matrix: 
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The plant system is easy to set up in our Mathematica 
package; a sample output is shown in Figure[3j The pack- 
age, based on the circuit modeling and analysis frame- 
work of Sarma et al. [32], allows one to arbitrarily con- 
catenate and link smaller elements to form larger quan- 
tum circuits, as long as all of the components are linear. 
The feedback control circuit is one example system the 
package can be used to simulate. 

Once the combined plant / controller system is set up, 
with its associated A, B, C and D matrices, the covari- 
ance matrix <xy = | (xiXj + XjXi) can be computed with 
the Lyapunov equation 



Act + crA T + BFB T = 
For a model system with the parameters 
k x = k 2 = k 3 = 0.01, A = 0.1 



(8) 



(9) 



we plot the cost function (a^a) as a function of noise k n 
for the various controller types in Figure |4j The orange 
line gives the performance of the plant without a con- 
troller. As expected, the photon number rises linearly 
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with the noise power. It is not hard to show that this 
matches the analytic result 



(ota) 
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that one can derive from the QSDE's. 

The trivial controller is simple enough that it also has 
an analytic solution. The two mirrors, rather than leak- 
ing photons separately, do so constructively so that the 
leakage amplitudes (rather than their powers) add up. 
This requires the replacement k\ + fc 2 — > (Vk~i + \fk~2) 1 
in (10), leading to the following result 



(11) 



fci + k 2 + k 3 + 2y/kik 2 
which agrees with the numerical data plotted in Fig. |4j 



A. Classical Controllers 

More sophisticated are the classical measurement con- 
trollers. The first simply makes a homodyne measure- 
ment of the ddix field. This signal is fed through a clas- 
sical circuit which generates an output. The heterodyne 
controller is slightly more complicated, and can be mod- 
eled as a two-input homodyne measurement controller in 
the following circuit (using the notation of j29j [33] ; see 
Appendix A): 



(Hom) 2 -m < (h ffl e l7r/2 ) < BS(a) 



(12) 
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FIG. 4: Bottom: Plant Ccavity photon number, as a function 
of noise strength fcjv. The uncontrolled case is shown, as 
well as the photon number for various control schemes. Top: 
Photon number relative to the no-control case. Smaller is 
better. 



In addition to the homodyne controller's parameters, we 
can also vary the beamsplitter transmittance. (Setting 
the beam-splitter transmission coefficient a — > 1 would 
send all the light entering controller input 1 into the x- 
quadrature homodyne detector, so the classical homo- 
dyne controller is really a special case of the classical 
heterodyne controller.) 

This example, in particular, illustrates the power of the 
Gough- James circuit algebra in treating control problems 
when the controller has a more complex, "circuit-like" 
structure. Having written code to output the ABCD 
model for a general n-input homodyne controller, it 
would have been straightforward, albeit tedious, to write 
additional code for the n-input heterodyne controller. 
But using the Gough- James circuit algebra allows us to 
write the n-input heterodyne system in terms of a 2n- 
input homodyne system, plus some beamsplitters and 
phase shifters, so we get the heterodyne controller for 
free. By breaking the system into smaller components, 
we can reduce the total amount of work we need to do in 
quantum control and simulation problems. 

There also exist "analytic" formulas for LQG-optimal 
classical controllers in the classical case. It is not diffi- 
cult to rewrite Eq. ^ in the standard form for an LQG 
problem: 



dx = Apxdt + Bpdu ■ 
dy — Cpxdt + dv 



dw 



(13) 



Here dy is the measurement signal, du is the controller 
output, and dw and dv are the plant and controller noises, 
dw ~ N(0,F w dt), dv ~ N(0,F v dt). Unfortunately, in 
this system the noises are correlated; the vacuum noise 
dcii acts on both the plant and, after reflection off mirror 
fci, the controller. One can define a covariance matrix 
= (dwidvk) to account for this correlation. 
A common trick is to remove the noise correlations by 
performing a change of variables [35] . Since dy — C p xdt — 
dv = 0, we can subtract this quantity from the first line 
of ( 13 ) to find an equivalent equation of motion: 



dx 



Apxdt + Bpdu + dw + MF V 1 (dy — C p xdt — dv) 
(A p - MF v C p )x + B p (du + B~ l MF~ l dy) 
+(dw- MF^dy) 



Ax 



Bpdu 



dw 



(14) 



Here, the noises dv and dw are uncorrelated. The con- 
troller for this plant will consist of a Kalman filter and a 
feedback: 

dx = (A- B p L - KC p )x + Kdy 

du = —Lxdt (15) 

The Kalman gain and feedback matrices can be obtained 
by solving the Riccati Equations: 

K = oC T F~ l (Aa + aA T - aC^F^Ca = 0) 
L = R- 1 B T \ (A T X + \A + Q-\BR- 1 B\ = 0) 

(16) 
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FIG. 5: Left: Optimal heterodyne amplification r\ as a func- 
tion of plant noise. Right: LQR as a function of controller 
amplification, for five different noise values. 



(Here Q and R are LQR optimization weights for the 
plant and controller states; we assume Q 3> R). For 
this particular case we optimized the classical controllers 
numerically, but the results agree with the analytical ex- 
pression. When optimizing the measurement controllers, 
we found that the best controllers always had dynamics 
that were much faster than the plant timescales. When 
this happens, the controller's internal dynamics can be 
adiabatically eliminated and the controller can be re- 
placed by a simplified "limit model" of the original com- 
ponent [3H 88J ■ When a linear component is adiabati- 
cally eliminated, its internal variables are removed and 
its ABCD model is replaced by the input-output rela- 
tions: 



da=(D- CA- 1 B)da 



(17) 



The homodyne controller, adiabatically eliminated, be- 
comes: 



da x = £ida x + daui, x 
da p = £ 2 da x + ddki, p 



(18) 



In this device, the signal da x is measured, amplified by 
factors ci and C2, and imprinted onto the output field. 
The downside of this measurement is the additional noise 
dcifci that the output accrues. 

The optimal heterodyne cont roll er uses a 50-50 beam- 
splitter so we set a = l/v2 in (12). It too has very fast 



dynamics that can be adiabatically eliminated to give: 



da x = £,(da x + da kltX ) + da k2)P 
da p = £{da p - da kl:P ) + da k2 , p 

Or cquivalcntly: 

dA = £{dA + dA\ x ) + dA k2 



(19) 



(20) 



The heterodyne controller amplifies both quadratures, 
but there is an additional noise due to splitting the beam 
before measurement, dA k i, as well as the measurement 
noise itself. The LQR can be computed analytically, and 
the analytic result agrees with the numerical optimizer. 
Setting £ = sinh(?7), we have: 



<at a > c 



k 2 sinh r\ + k^k n 



k\ + k 2 



2y/kik~2 sinh.77 



(21) 
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FIG. 6: Optical parametric oscillators adiabatically eliminate 
into ideal squeezers. 



This is plotted in Fig. [5j As the plant noise increases, so 
does the controller's optimal amplification. It does not do 
well to increase the amplification indefinitely, however, 
since this also adds noise into the system. From Fig. 
[4j one can also see that measurement control does well 
at reducing the photon number for large k n , but in the 
quantum regime, k n < 1, it has hardly any effect at all. 



B. Coherent Control 

The three coherent controllers of interest are the cavity 
controller and the two OPO setups, as shown in Figure]!] 
The optimizer consistently showed that the best cavity 
controller is in fact the trivial controller (which is the 
special case of a cavity with mirror transmittivity set to 
zero). Because of this, we do not consider empty cavity 
controllers in this section. The OPO controllers, on the 
other hand, have more interesting behavior. 

As in the classical case, it was discovered that the best 
coherent controllers always had dynamics that were much 
faster than the plant timescales and could be adiabati- 
cally eliminated. A single OPO will adiabatically elimi- 
nate to a squeezer with the following input-output rela- 
tions: 



da T — e v da 7 



da r . 



~ v da r . 



(22) 



(up to input and output phase shifts). An OPO system 
with squeezed inputs and outputs, which can in princi- 
ple replicate any 2-port linear quantum system with a 
single internal degree of freedom [33], will adiabatically 
eliminate to arbitrary two-mode squeezer (Fig. [6]). As far 
as this control problem is concerned, the best two-mode 
squeezer is the linear amplifier, given by the input-output 
relations: 

ddi = cosh(?y)(iai + sinh(r/)c?a2 

da 2 — smh(rj)da 2 + cosh(r])da 2 (23) 
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Analytic formulas can be derived straightforwardly 
from the quantum stochastic differential equations. For 
the squeezer: 



(ata> s 
where 



Re 



(k 2 sinh 2 T] + 2k n ) - 2 



k-2\/k\k2 cosh 7/ sinh rj 
G+2iA e 



Re [G - ±k x k 2 sinh 2 rj/(G + 2iA)] 



(24) 



G = ki + k 2 + /J3 + 2\J k\k 2 cosh(?7)e 



For the linear amplifier: 



(a)a) 



k 2 sinh 77 + fc3/c r , 



2-sq 



2\/k\k 2 cosh 77 



(25) 



Qualitatively, the results for the heterodyne controller, 
Eq. (21 ) and the linear amplifier, Eq. ( [25] ) look very sim- 
ilar. Both the heterodyne controller and the linear am- 
plifier reduce the cavity's photon number by amplifying 
the feedback signal, but also add noise to the system. For 
equivalent levels of amplification (compare (21 1, substi- 
tuting sinh ry —> cosh?/, to (25)) the classical controller 



adds extra noise into the system from the measurement 
process. When k n and 77 are large, this extra noise is 
negligible, but in the quantum regime where k n and 77 
are < 1, this noise can play a major role in making the 
linear amplifier outperform the heterodyne controller. 



As far as optimization is concerned, Equations (24-25) 
are simple enough to apply. Finding the best controller 
just involves minimizing these functions with respect to 
77. But remember that it was not at all obvious that 
the best quantum controller should be an adiabatically 
eliminated squeezer. This had to be demonstrated by 
optimizing the general OPO controller, which has many 
more parameters, and comparing the result to that of 
the squeezer. This required a Mathematica package to 
quickly convert circuit diagrams to ABCD models, and 
an efficient optimizer to find the best controller parame- 
ters. 

Notice that, for large k n , the performance of the two 
quantum controllers follows the classical performance. In 
the classical limit, the OPO / squeezer is amplifying a 
single quadrature and feeding this back into the plant 
(with the proper phase shift). Likewise, the classical 
controller measures a single quadrature, amplfies that 
signal and sends this back into the plant. Thus, the 
OPO / squeezer is a "homodyne-like" controller in the 
classical limit. By contrast, the linear amplifier amplifies 
both modes equally and feeds back the result, making 
it a "heterodyne-like" controller which tracks the perfor- 
mance of the heterodyne controller in the classical limit. 

However, in the quantum regime, this distinction is 
lost and both of the quantum controllers track the perfor- 
mance of the trivial controller. Below a threshold value 
of 



2^k 2 



\fkik 2 kz 



(26) 
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FIG. 7: Left: Optimal squeezing for the squeezer (solid) and 
two-mode squeezer (dotted) controllers. Right: Performance 
as a function of squeezing for multiple noise levels. 



(for this system, k n mm — 5), any squeezing will increase 
the noise in the cavity, so the optimal value of 77 is zero 

in other words, for k n < 5, the best controller is the 
trivial controller. 

As Fig. [7] illustrates, when k n > 5, the best controller 
has a nonzero amount of squeezing. We plot the con- 
troller performance as a function of squeezing for five 
different noise levels on the right pane of the figure. In- 
tuitively, this is a battle between the noise introduced 
by squeezing and the noise removed by constructive in- 
terference with the light leaking out mirror 2. When r\ 
is low, the latter dominates. By increasing the squeez- 
ing, we effectively increase the amplitude of the field 
impinging upon mirror 2. Recall that the trivial con- 
troller worked by constructive interference between this 
field and the light leaking out of mirror 2. By increas- 
ing this field's amplitude, we magnify the effect of this 
interference; this reduces the overall cavity photon num- 
ber. This explains the cosh r\ term in the denominator of 
( p5| ) . But a squeezed vacuum carries photons of its own, 
and some of these photons leak back into the cavity. If 
the squeezing is too high, this winds up increasing th e 
photon number, giving rise to the sinh 2 77 term in (25). 
Above the threshold temperature fc n ,min, the ideal 77 lies 
somewhere between these extremes. 



Below the threshold temperature, the cavity photon 
number is so low that the interference effect never wins 
out - squeezing the control field always introduces more 
photons in the cavity, and the best controller involves 
looping the output from mirror 1 into mirror 2 without 
squeezing the trivial controller. 



III. OPTICAL FEEDBACK CONTROL OF A 
MECHANICAL OSCILLATOR 

Optomechanical oscillators - mechanical springs that 
couple to an optical field via a cavity - have been a 
topic of tremendous recent interest in the physics com- 
munity [39]. A central goal has been to find ways to 
exploit optomechanical coupling to cool the mechanical 
oscillator from ambient temperature to its ground state, 
using optical feedback. 
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In this section we analyze the optomehcanical oscilla- 
tor as a coherent control system, with the spring compris- 
ing the plant, and with optical probing and feedback. We 
optimally cool the oscillator by solving the LQG con- 
trol problem for the cost function (b^b), where b is the 
spring's annihilation operator. While the control setups 
we consider may appear impractical from an experimen- 
tal perspective, we will discuss how they can be related 
to systems that are more realistic to implement. 

At the heart of this control problem is the "adiabati- 
cally eliminated cavity," depicted in Figure [HJ If we go 
into the rotating frame for the light, this has the SLH 
model 



S = 1 



2x2, 



L 



fna, y/fl/Qb 



h = hntfb- 



-•qa' ax m 
(27) 

where 572 is the natural spring frequency, Q is the Q-factor, 
k is the cavity decay parameter, and m is the mirror 
mass. See Table [H 

System ( |27[ ) is nonlinear by virtue of the interaction 
term rja^aX. This term is due to the photon pressure of 
the field in a cavity, which exerts a physical force on the 
mirror. In the limit that the light mode a evolves much 
faster than the mechanical mode b, we can adiabatically 
eliminate the former to give an SLH system of the form: 



S = 
where 



e i<t>{x m —x m0 ;r]/K) q 

1 



, L = 



{z-.rj/n) = 2tan (2t]z/k) 



(29) 



is the phase shift of the cavity reflected light, as a func- 
tion of the mirror position (we have absorbed a factor —1 
in S for convenience) . This is still a highly nonlinear sys- 
tem. A real optomechanical oscillator is usually driven 
by a coherent field, and the output that is measured is 
generally interfered with an equal and opposite field, so 
as to discern the phase fluctuations on a homodyne de- 
tector. Thus, the real plant system we are interested in is 
the adiabatically eliminated cavity sandwiched between 
two coherent displacements. For a cavity subject to a 
coherent input of amplitude r±, we write this as: 

(Cavi) = L(-r) < (Cav) < L(r) 



This has the simple, linear SLH model: 

Si = l 2 x2, Li = [K lXm , v^VQfe] , Hi = Qtfb (31) 

with Ki = Arjri/ Ki is the effective coupling between the 
spring and the field, which need not be positive or even 
real. The x m -coupling to the field da\ gives rise to the 
following input-output relations: 



(dx m )i 
(dp m )i 







— 2K\dai v 

ddi x = dcii x + 2KiX m dt 
do,\p — dd\p 



(32) 



0, y/tt/Qb , H = ntfb 
(28) 



The state variable x m is imprinted on the output da\ x , so 
by measuring the x-quadrature of the output field, we can 
deduce the value of x m ; this allows us to use the mirror as 
a "measurement" device, learning information from the 
output field. Note that this only works for ddi x ; no infor- 
mation is imprinted onto the p-quadrature of the output. 
Conversely, by sending in a particular input da\ p , we can 
alter the state of the system; this allows us to use the 
mirror as a "feedback" device. Note likewise that feed- 
back is not possible via the da± x channel, which does not 
affect the system. 



A. Plant System 

The plant-controller setup is shown in Figure [9] The 
plant system consists of two (adiabatically eliminated) 
cavities coupled to the same mirror. The output from 



Qty 


Value 


Ki 


= Arjn/Ki 


n 


= y/Pi/hu 


K i 


= Uc/2U 


V 






= (1 - e -«V fcT )-i 


km 


= fi/Q 




(30) 


Qty 


Description 


Typical Values 




P 


Laser power in coherent displace- 
ment n, i — 1,2 


1 pW-1 mW 


dA 2 

.r 


U 


Power transmittance for cavity mir- 


io- 6 -io -3 


k 


ror i. Inversely proportional to fi- 
nesse. 

Length of cavity i 


10" 6 -10 _1 m 




m 


Mass of spring-mounted mirror 


Kn 15 -l(T 10 kg 


"dA 2 


n 


Spring oscillation frequency 


kHz-GHz 


Q 


Spring quality factor 


10 3 -10 7 






Laser frequency 


2-4 x 10 15 /s 



FIG. 8: Single cavity with modes o, a 1 , coupled to a mechan- 
ical oscillator with modes b, b. 



TABLE I: Parameters for the optical cavity controller prob- 
lem. See, e.g. [3S] 



x m , p m: or b, b* 




Measurement Simple Cavity OPO 
Controller 



FIG. 9: Control-system setup for the mechanical oscillator 
cooling problem. Four potential controller designs. 



dA^ — |Cav 2 
dA> 



Cav! 



dA, 



Spr 



-dA x 

dA k 

-dA, 



FIG. 10: Equivalent view of the plant-controller setup shown 
in Figure [9] See Eq. (33). 



the first cavity, dA\, goes into the controller, and the 
controller output is fed back into the second cavity in- 
put oL4 2 . Not shown are the two coherent displacements 
(lasers) putting fields dA\ and dAk into nonvacuum co- 
herent states. These coherent fields allow us to replace 



the cavity with model ( 28 1 with the linearized model ( 31 1 . 



Since the system is now linear, this becomes an LQG con- 
trol problem. The combined plant-controller system can 
be viewed as a feedback loop from output dA\ to input 
dA2, or conversely, we can write it as a series product 

(Sys) = [(Cav 2 mh)<K< (Cavi ffl h)] ffl (Spr) (33) 

where (Sys) is the combined system, Cavj is the i th cav- 
ity, with SLH model (1, \/kiX m , _), (Spr) gives the spring 
and phonon couplings, SLH model (1, \/k m x m , flwb), 
and AC is the controller. See Figure [T0| 

The controllers we consider here are not unlike those 
for the simple cavity. It is not difficult to show using 
Eq. ( 33 ) that the trivial controller amounts to no con- 
trol at all at best, and additional noise at worst. The 
classical controller measures the output from mirror cav- 
ity 1 and sends an input in to cavity 2, as a function 
of the controller's internal state. (Note that we only 
need to consider a classical controller that measures the 




No Controller 
Classical Homodyne 
Homodyne, K t » K ; 
Simple Cavity 
OPO Controller 



10 5 




FIG. 11: Bottom: Plot of the average phonon number (N) = 
(6 6) of the mechanical oscillator for three different control 
schemes. Top: Phonon number reduction, relative to no- 
control case. The general coherent controller result is not 
shown, since it overlaps the OPO line, the optimal coherent 
controller being an OPO cavity. 



x quadrature dd\ x ] dd\ p contains no information about 
the plant's state.) The simple cavity and OPO cavity 
coherently process the signal rather than destroying it in 
a measurement. Finally, we considered the most general 
coherent controller, an open quantum system specified by 
arbitrary A,B,C,D matrices satisfying the realizability 
relations. For the LQG problem of minimizing (b^b'), we 
found optimal controllers in each class for the following 
plant system: 



SI = 100 (arbitrary units) 

k rn = 0.01 

Q = 10000 

k„ = 10~ 9 -10 9 



(34) 



In the optimization, we are allowed to vary both the 
controller parameters and the couplings K\, K% to the 
cavities in plj ). This is because the couplings depend 
on the input laser powers Pi (in addition to the mir- 
ror transmittances ti), which are external quantities (see 
Table |l| rather than fixed properties of the plant itself. 
Here we will operate primarily under the assumption 
K\ = —Ki = K; this is a reasonable assumption that 
avoids classical solutions with divergent controller gain, 
but we also show that the coherent controllers discussed 
here outperform the best classical controllers even when 
this assumption is relaxed. 



Figure 11 plots the performance of the measurement, 
simple cavity, and OPO controllers. For very low ambi- 
ent temperatures where the noise is weak, the plant is 
nearly in its ground state to begin with, and none of the 
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controllers can reduce its value. This differs from the op- 
tical cavity. In the cavity, we used a "trivial controller" 
to cause the light leaking out of mirror 1 to interfere con- 
structively with the light leaking out mirror 2, increasing 
the net dissipation from k\ + k 2 to (\fk\ + v^) 2 - No 
such scheme exists in the oscillator because phonons do 
not "leak out" of the system the same way photons leak 
out of an optical cavity. 

At high temperatures, the best classical controller 
and the OPO controller do equally well, each reduc- 
ing the phonon number by a factor of exactly Q = 
10000. The cavity controller does reasonably well, reduc- 
ing the phonon number by a factor of about 0.354Q = 
3540. These results are not very surprising. The high- 
temperature limit takes our oscillator into the classical 
regime, where vacuum noise is negligible and no coher- 
ent controller can hope to outperform the best classical 
controller. 

The interesting region lies between these two limits. 
Here, there is a sharp cutoff, near k n 0.2, below which 
the classical measurement controller becomes useless. As 
explained below, the classical controller must add noise 
to the system to make a measurement; below a certain 
threshold, the gains from control are offset by the noise 
from measurement. In this region, the cavity and OPO 
controllers do significantly better than the classical con- 
troller, in some places by a factor of 100-200. 



B. The Classical Controller 

The classical controller works by measuring the plant 
output field (dA in Figure 12) and inferring the plant's 
state from this measurement. From the inferred plant 
state, the controller applies a feedback signal, which is 
added to an auxiliary vacuum input dAk and sent back 
to the plant. 

The plant output contains two quadratures, but only 
one of them contains information about the system. 
Thus, in our classical controller we choose to measure 
the ^-quadrature of the output, and necessarily discard 
the dA p . This is the optimal control strategy in the clas- 
sical case because dA p does not contain any information 
about the system. Like any LQG-optimal controller, the 
classical controller consists of a Kalman filter, which es- 
timates the plant state, plus a feedback element. 




Coherent 




The classical controller adds two sources of noise to the 
plant. First, by sending a laser through the measurement 
cavity (Cavi), it adds measurement noise, with an am- 
plitude that scales as O(K). Second, the feedback field 
dA (with a vacuum noise component due to the auxiliary 
field dAk) is sent through the controller, adding a feed- 
back noise of equal magnitude, also O(K). Both of these 
factors increase the cavity phonon number by 0(K 2 ), in- 
dependent of the noise k n . The control loop will decrease 
the cavity phonon number by an amount proportional to 
the present phonon number, which increases with k n . In 
the high-fc„ limit, the "control" term dominates and the 
coupling K is large. By contrast, in the low-fc„ limit, the 
"noise" term is dominant, and the optimal value of K is 
small or zero - no measurement controller can effectively 
reduce the phonon number, since the noise incurred will 
more than offset any gains from control. 

An important thing to note is the role the p-quadrature 
field da p plays in this noise budget. It is true that da p 
does not contain any information about the plant state. 
But this quadrature still plays an important part, since 
da p gives rise to the noise in the measurement cavity, 
and dd p gives rise to the noise in the feedback cavity. 
Because da p and dd p are independent (the former being 
destroyed in the da x measurement), their noises add up. 
The beauty of coherent control is that we can process the 
x-field without destroying da p and the measurement and 
feedback noises become correlated. If done right, they 
cancel each other out. 

If we are free to relax the K\ = —K% assumption, then 
the classical controller does somewhat better (dashed 
line in Figure 11), but still undcrpcrforms the coherent 
schemes discussed below. When K\ ^ K 2 , the optimal 
classical controller tends to have K 2 <C Ki , which greatly 
suppresses the measurement noise. To compensate for 
this disparity, the controller must have a large classical 
gain. 

It might be thought that a heterodyne-based control 
scheme, like that in Figure [13] could perform better than 
the best homodyne controller. After all, the homodyne 



dA 



k2 



du dA 7 dA, 




1 — — \ 



Classical 



FIG. 12: Flow of x- and p-quadrature signals (blue and red, 
respectively) in the classical and coherent controllers. 



FIG. 13: Heterodyne-based measurement controllers, which 
measure both quadratures of the beam by splitting it, do not 
not outperform the best homodyne controller for this system. 
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controller is just a special case of the heterodyne con- 
troller, where the beamsplitter has a transmissivity of 
100%. Moreover, one might imagine using a heterodyne 
scheme to cycle part of the da p quadrature back into the 
plant, canceling out part of the measurement noise with 
the feedback noise. However, we find numerically that 
the most general heterodyne controller does not perform 
any better - either with K\ = Ki or not. The extra 
noise added from splitting the beam outweighs any of 
the benefits of the control scheme. 



C. Simple Cavity Controller 

An empty optical cavity with two input / output ports 
has the following SLH model: 

S = 12x2, L = [y^Kia, \/K2a] , H — Aa'a (35) 

Here the re's are mirror decay parameters and A is the 
detuning of the cavity. The QSDEs for the cavity are 
easy to derive: 



da = (— iA — n/2)a dt + ^/kicLAi 
dAi = dAi + J~K\a dt 



K 2 dA 2 



(36) 



Remember that, in addition to the controller parameters, 
we can vary the input coherent fields, which allows us to 
vary the plant's a; m -coupling K . The laser field impinging 
on cavity 1 adds shot noise to the mirror; in this setup, 
since K\ = —K2 = K, the shot noise from cavity 1 will 
exactly cancel the shot noise from cavity 2 (if we let K\ 
and K 2 vary freely, we find that the optimal controller 
has K\ = — K2). As a consequence, the cavity controller 
has neither measurement nor feedback noise. 

The optimal detuning and couplings are plotted in 



Figure 14 Not surprisingly, as the noise on the mir- 
ror is increased, the couplings if 1 2 and k 12 increase as 
well. The detuning A, which shows no dependence on 
the noise power, always remaining at a constant value 



A « il — 100 for this system, making the cavity con- 
troller setup analogous to two coupled harmonic oscilla- 
tors, one mechanical and the other optical [JO]. Absent 
the couplings, the quadratures x = a + ,p = (a — <v)/i 
would evolve just as the mirror variables x m ,p m . 

This can also be interpreted as a form of sideband 
cooling. The detuning A « indicates that our con- 
trol system is being driven by laser light at a frequency 
w cav — fi, where ui cav is the cavity resonance frequency. 
The plant-controller coupling serves to convert photons 
of frequency w cav — Q to photons of frequency u; cav , cool- 
ing the oscillator. At high temperatures, we need a large 
cooling rate to counter the noise; this is achieved by using 
a cavity with a broad bandwidth re, so that both oj cav — f2 
and w cav photons are interact effectively with the cavity. 
Conversely, at low temperatures, we need to work in the 
resolved sideband limit re <C to suppress quantum fluc- 
tuations of the radiation-pressure force [4"3Tt45| . 

The effects of this cooling are made manifest on the 
output power spectrum of the photon channel Pi(u>) — 
Ai(uiy Ai(u>), where A\{uS) is the Fourier transform of 
the stochastic process dA\(t). In the frequency domain, 
the relevant QSDEs for the combined plant-cavity system 
are 



(-iA - (ki + re 2 )/2)a + ^/KiK(b + b r ) 
+iojy/KiAi + iwy / re2A 2 



[(-ifi - n/2Q) b - y/^K(a - a f )] + iujy / n/QA 3 
—iuiAi + -y/rela (37) 



—iuib 
-iwAx 



This power spectrum is plotted in Figure |15| As the 
exiting light is blue-detuned, it reduces the phonon num- 
ber in the oscillator, driving it towards the ground state. 
For small k n , when the plant and controller are weakly 
coupled, there is a single sideband corresponding to the 
plant's oscillation frequency When k n is large, the 
plant and controller become strongly coupled and the 
combined system resonates at two different frequencies, 



A (solid) 



10 - 




A - ,, AT, (bold) 



k 2 (dolled) 



/ / K\ (dashed) 

0.001 - I 



10- 5 L 

™j ^ _J _J _J l n 

io-» 0.001 1 1000 10' 10' 

FIG. 14: Parameters of the optimal simple cavity controller, 
as a function of noise strength. 




FIG. 15: Calculated output spectrum of light exiting the op- 
timal simple cavity controller. Six values of k n are plotted, 
10 5 (darkest), 10 3 , 10 1 , HT 1 , 10" 3 , and 10" 5 (lightest). 
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one larger than f2 and one smaller. This is the origin of 
the sideband splitting in the figure. 

The system can also be understood as a form of coher- 
ent Kalman filtering. Recall that the optimal classical 
controller works as a Kalman filter, reproducing the state 
of the plant by measuring one of its outputs. The cost we 
paid for the Kalman filtering was additional noise added 
to the system. The cavity controller can also be thought 
of as a Kalman filter, but one that preserves the coher- 
ence of the input signal dA. From a quantum mechanical 
standpoint, in the classical controller, the p-quadrature 
da p is essentially discarded after the measurement. In the 
cavity controller, the field retains its coherent properties 
and the da p coming out is the same as da p going in. This 
makes the noises in the measurement and feedback cavi- 
ties correlated. In the present setup, they exactly cancel 
out. This cancellation of the measurement noise is what 
gives the coherent cavity controller its superior perfor- 
mance, particularly in the low phonon-number regime. 

Measurement sensing experiments |41j . particularly in 
the context of LIGO [42], show similar improvements, 
but for a different performance metric. This suggests 
that LQG control is far from the only problem to ben- 
efit from this noise cancellation and coherent feedback; 
similar gains should be expected in all types of control 
problems when the plant operates in the quantum regime. 



D. OPO Cavity Controller 

Recall from Eq. ^ that the OPO has the following 
SLH model: 



S = 1 2 X2, L = [y/Kta, y/j^a] 



H=-x T 
4 



Aa f a + 



A - Im(e) Re(e) 
Re(e) A + Im(e) 

e*a 2 -e(a)) 2 
2i 



(38) 




10-» 10-" 0.001 1 1000 10" 10' 

FIG. 16: Parameters of the optimal OPO cavity controller, 
as a function of noise strength. 



For fullest generality, the OPO controller is placed be- 
tween two phase shifters, so the actual controller is 
e i<t>\ <] (OPO) < e^ 2 . Between the controller, the phase 
shifters and the couplings K% t 2, there are nine free pa- 
rameters in this LQG problem. The best OPO controller 
parameters, found using the optimization code, are plot- 
ted in Figure [16] As with the cavity controller, the best 
OPO controller has K\ = -K 2 . 

For k n < 1800, the OPO behaves much like the simple 
cavity. Its detuning is close to f2, the coupling K\ — 
—K2 increases with k n , and the mirror losses K\, k%, while 
small, increase with increasing noise (k 2 is too small to 
be seen on this plot). For the most part, e« A and the 
OPO squeezing is only a perturbation on the dynamics 
of an empty cavity. 

At k n « 1800, this changes suddenly. This happens 
because the OPO controller has two local minima. Be- 
low k n ~ 1800, the empty cavity-like local minimum is 
smaller, but above this threshold, a new minimum dom- 
inates. In this regime, the coupling K is much stronger 
than before and the mirrors k± , k 2 are much more lossy. 

The OPO controller appears to be the best coherent 
controller one can make for this system. We ran the op- 
timizer for general coherent controller, subject to no con- 
straints other than the realizability conditions Q . At no 
point did we find a coherent controller that outperformed 
the OPO for this system. This in mind, ths discontinuity 
at k n can be better understood. As the best relizable 
controller, the OPO must do at least as well as both the 
simple cavity and the classical controller. For weak noise, 
the simple cavity outperforms the classical controller, so 
we expect the OPO to look more like a simple cavity. 
For strong noise, the classical controller does better, so 
we expect the OPO to look more like a classical con- 
troller, inasmuch as this is possible. There is no reason 
to assume that the transition between the two must be 
smooth. It may be marked with bifurcation points, as in 
Figure [7] for the cavity control problem, or it may occur 
with a discontinuity in the parameters. What happens 
for a general plant / controller system will depend on 
the landscape of the cost function, and in particular, the 
behavior of local minima. 



E. More Realistic Control Systems 

The control systems discussed above can be imple- 
mented in principle, but they require two separate mir- 
rors and two separate cavities to be coupled to the same 
mechanical oscillator, which may prove difficult to build 
in a laboratory. Fortunately, one can show that for the 
cavity controller and the OPO controller, equivalent sys- 
tems can be realized using a non-adiabatically eliminated 
cavity with one of its mirrors on a spring. 

First, the simple cavity controller. Recall from (33) 
that the cavity controller system can be modeled as 



[(Cav 2 ffl h) < (Cav) < (Cavi ffl h)} ffl (Spr) (39) 
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^nv Prr 



dA, 



dA, 



FIG. 17: Model for a non-adiabatically eliminated cavity. 




FIG. 18: Model for a non-adiabatically eliminated OPO cav- 
ity with a spring mirror. 



which has the SLH model 

S = 12x2, L = [yfj^a, y/K^a, y/n^b] 

H = u} c a^ a + Qb^b + y /r KiKix m p c (40) 

Now consider a system, depicted in Figure |17[ consist- 
ing of a non-adiabatically eliminated cavity with one of 
its mirrors attached to a spring. This has the SLH model: 



5 = 13x3) L = jKia, Jn^a, V k m b 



H = A (Ja + ntfb 



77a' ax m 



(41) 



A laser L(a) sends a coherent input into mirror 1, giving 
the system Cav < (L(a) ffl 12)- Of course, the internal 
dynamics do not depend on anything downstream of the 
system, so we can just as well use (L(a') ffl I2) < Cav < 
(L(a) EE I2), for any a' . Making substitutions a — > a — 
ao, b — > b — bo to center around the equilibrium point, the 
SLH model becomes: 



S — 1-3x3) L — 

H = Aa f a + fltfb 



K\a, Y / K2~a, y k m b 



A 2 + (k/2) 



2 '$-m*&c 



(42) 



Ignoring the nonlinear term, this is almost identical to 
(40 1. One can convert the x m x c term to an x m p c term 



with a canonical transformation, and the coefficients can 




Vacuum 



K T vac = 


System 

l~sys 
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^ K T 5ys 


k m T h 



k m = O/Q 



FIG. 19: Coherent control problem represented as two cou- 
pled thermodynamic systems. 



be matched by varying a. Thus the systems in (40) and 



(42 1 arc equivalent, and the "simple cavity controller" 
can be realized in the lab using a single cavity with a 
mirror attached to a spring. Cooling an oscillator in this 
setup has been realized experimentally, though it was not 
interpreted as a control system [TTHT5] . 

The OPO controller is just like the cavity controller, 
but the Hamiltonian has an additional squeezing term; 
see ((5j). The same procedure can be applied to show 
that the OPO plant-controller system is equivalent to a 
(non-adiabatically eliminated) OPO cavity with a spring 
mirror, as shown in Figure [T8| 



F. Quantum Refrigerator Analogy 

One thing we notice from the optimal controller per- 
formance is that, in the strong-noise limit, the optimal 
controllers - classical, OPO, cavity - all reduce the spring 
phonon number by a factor of about Q = Q/k m . The 
classical and OPO controllers reduce it by exactly Q, 
while the cavity controller only reduces it by a factor 
Q/2.83. This factor-of-Q reduction can be understood by 
viewing the plant and the controller as thermodynamic 
systems. 

Figure [T9] illustrates our point. Starting with a cavity 
with a spring mirror, we separate the system into the 
cavity, which oscillates at a frequency w c , and the spring, 
which oscillates at a frequency ft. Each system has its 
own coupling to the environment. The cavity couples 
to a vacuum-state environment (T — 0) with coupling 
strengths K\, K2, the spring, couples to a heat bath with 
T/j > with strength k m , and a spring-cavity coupling 
K\ = —K 2 couples the two modes. 

If the spring and cavity oscillate at about the same 
frequency and the spring-cavity coupling is strong com- 
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pared to the other two, then the "temperature" of the 
spring will be roughly equal to the "temperature" of the 
cavity. We denote this temperature T sys . One expects 
the combined system to be in thermal steady-state with 
both the heat bath and vacuum inputs and outputs; this 
gives us the energy balance equation: 



(43) 



where k = K\ + k 2 ~ u c ~ SI, and k m = fl/Q. Solving 
for the system's steady-state temperature, 



sys 



(44) 



From general arguments, we can therefore expect that 
most good controllers will reduce the spring phonon num- 
ber by a factor of about Q, but that no controller will 
do significantly better. Note that, since this argument is 
based on thermodynamic assumptions that are only ap- 
proximately valid here, the factor-of-Q reduction is only 
approximate, and only holds in ths classical limit. These 
classical results, unsurprisingly, break down in the quan- 
tum regime because, among other things, the effects of 
vacuum noise inputs become important. 



IV. CONCLUSIONS 
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FIG. 20: The concatenation, series, and feedback products are 
the three basic operations of the Gough- James circuit algebra. 



in linear systems may provide important clues for the de- 
sign of quantum controllers for nonlinear systems such as 
optical switches or error correcting codes. 
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In this paper, we have studied the coherent-feedback 
cooling of linear quantum systems from an LQG con- 
trol perspective. The systems were modeled using the 
SLH framework and the Gough- James circuit algebra, 
which allow arbitrarily large circuits be constructed in 
a straightforward and systematic manner. The evolu- 
tion of the system was studied using QSDEs, the open- 
system analogue to the Heisenberg Equations. We wrote 
Mathematica scripts based on the QHDL/M framework 
to model quantum LQG control systems, and designed al- 
gorithms to optimize a controller's parameters for a given 
setup. 

For any LQG control problem, there is always a quan- 
tum controller that does at least as well as the optimal 
classical controller. In the quantum regime, when ex- 
citation number in the plant is of order unity, we have 
shown that the best quantum controller can do better - 
in some cases, significantly so. Two systems - the optical 
cavity and the optomechanical oscillator - were studied 
in detail. For the former, modest gains were found using 
coherent control in the low-photon-number regime. For 
the latter, the gains were much larger. 

One could imagine extending these results to look at 
non-quadratic cost functions in linear control systems. 
Indeed, some work has already been done on this mat- 
ter, focusing on using coherent feedback to maximize the 
squeezing in a cavity mode [17] . Taking a control theory 
perspective may also provide insight into minimizing the 
noise in optomechanical sensors. In addition, the under- 
standing the superior performance of coherent feedback 



Appendix A: SLH, QSDE, ABCD models 

An open quantum system with n bosonic, Markovian 
input-output channels, can be represented as a triple [28, 

(S,L,H) (Al) 

where S is a unitary n x n operator-valued matrix, L is 
an n-component operator-valued vector, and H is Her- 
mitian. 

A quantum circuit is built up by connecting together 
components of this form. Any such a circuit can be writ- 
ten in terms of a circuit algebra of concatenation, series 
and feedback products, shown in Figure [20| 

The concatenation product joins two systems with- 
out connecting any of the input/output ports. If G\ — 
(Si,Li,Hi) and G 2 = (S 2 , L 2 , H 2 ) then the concatena- 
tion is given by: 



Si 




Li 










s 2 




_L 2 _ 


, Hi + H 2 j 



The series product feeds ths outputs of one system into 
the inputs of another, and has the following SLH model: 

G 2 < Gi = (s 2 Si, L 2 + S 2 L X , H 1 +H 2 + Im(4s 2 Li)) 

(A3) 

The feedback product can [G]i->.j corresponds to taking 
output i of G and feeding it back into input j. This can 
also be represented with an SLH model. See Ref. 33 . 
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The Quantum Stochastic Differential Equations (QS- 
DEs) are Heisenberg-picture equations of motion for open 
quantum systems. They relate the evolution of the inter- 
nal state variables (denoted X) and the output fields dAi 
and gauge processes dAy to the inputs dAi, dAij, where 
the inputs are vacuum quantum Wiener processes. For a 
given SLH model, the equations are |29j : 



dX 



-i[X,H] + - 

+dA\sl j [X,L j ] 
[L],X]S ji dA i + 



L][X, Li\ + [L],X]L, 



dt 



S ik XS kj ~ X*ij ) dAij (A4) 



dA, 
dAa 



SijdAj 



Ljdt 



(A5) 



S* k dA kl S% + S* k dAlL., + L]dA k Sl + L\L 



(A6) 



Likewise, the Master Equation is the Schrodinger-picture 
equation of motion for an open quantum system. It gives 
the evolution of the system's density operator: 



f - i[p, H] + U lP L\ - \h\Up - \pL\u 



(A7) 



We form matrices S and A, and vector A by stacking S, 
A, and A: 



where 



Snh — 



A q = 



2M f 
2A/ f 
2M^ 



Sab 

o s* ab 



A a 

a: 

Aa 

A* 



Jlny.1n — In® 



1 

-1 



(A10) 



(All) 



is the canonical antisymmetric matrix of dimension 2n 
(written above as J, where the dimension is inferred), 
and 



A linear system has a quadratic Hamiltonian and lin- 
ear environment couplings. It takes the following SLH 
model: 



e. . — c. . 

^>ij — ^ij 

Li = AjX + A.; 

H = -x T Rx + r T x 



From these, the ABCD matrices take the form: 

A = 20 (r 
C = A, 
a = 26 (r 



jA T JA 



|A T JA 



B = 6A T JS, 
D = S, 
c = A 



(A8) 



(A9) 



M2nx2n — In 



(A12) 



The matrix S is made from the blocks S a b above, and 
likewise for A and A. 

These are similar to the formulas used in [TB], the dif- 
ference being that our input and output fields da are 
written in the quadrature basis, (da x ,da p ) rather than 
(da, da)). This keeps all of our inputs and outputs Her- 
mitian, and also makes all of our matrices real. 
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